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We study the topological charge and the topological susceptibility in lattice QCD with two de- 
generate flavors of naive Wilson fermions at two values of lattice spacings and different volumes, 
for a range of quark masses. Configurations are generated with DDHMC/HMC algorithms and 
smoothened with HYP smearing. We present integrated autocorrelation time for both topologi- 
cal charge and topological susceptibility at the two lattice spacing values studied. The spanning 
of different topological sectors as a function of the hopping parameter K is presented. The ex- 
pected chiral behaviour of the topological susceptibility (including finite volume dependence) is 
observed. 
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Figure 1: (a) Topological charge versus smearing levels for APE and HYP smearing at K" = 0. 15825, j3 = 5.6 
and volume 24^ x 48. (b) Topological susceptibility versus HYP smearing levels at JC = 0.15825, /3 = 5.6 
and volume 24^ x 48. 

1. Introduction 

Earlier attempt [1] in lattice QCD to verify the suppression of topological susceptibility with 
decreasing quark mass (m^) with unimproved Wilson fermions and HMC algorithm was unable to 
unabiguously confirm the suppression. Since topological susceptibility is a measure of the spanning 
of different topological sectors of QCD vacuum, the inability to reproduce the predicted suppres- 
sion may raise concerns about the simulation algorithm and the particular fermion formulation to 
span the configuration space correctly and/or efficiently. In this work, we perform a systematic 
study using unimproved Wilson fermions and demonstrate the suppression of topological suscepti- 
bility with decreasing quark mass. This work is part of an ongoing program [2, 3] to study the chiral 
properties of Wilson lattice QCD. For details on measurements of topological charge and suscepti- 
bility, and numerical values of measured quark masses, pion masses and topological susceptibilities 
see Ref . [4] . Our results of the topological susceptibility favourably compare in Ref . [4] with that 
from a recent mixed action calculation employing clover fermions and overlap fermions for the 
sea and the valence sectors respectively [5]. The detailed account of low lying spectroscopy will 
appear separately. 

2. Measurements 

We have generated ensembles of gauge configurations by means of HMC [6] and DDHMC 
[7] algorithms using unimproved Wilson fermion and gauge actions with «/ = 2 mass degenerate 
quark flavours at two values of the gauge coupling (j3 =5.6, 5.8) . At /3 = 5.6 the lattice volumes ai^e 
16^ X 32, 24-^ X 48 and 32^ x 64 and the renormalized quark mass varies from 15 to 100 MeV (MS 
scheme at 2 GeV). Quark masses are determined using axial Ward identity. At j3 = 5.8 the lattice 
volume is 32^ x 64 and the renormalized quark mass ranges from 20 to 90 MeV. The lattice spacings 
determined using Sommer pai^ameter at j3 = 5.6 and 5.8 are 0.077 and 0.061 fm respectively. 
All the configurations for lattice volumes 24^ x 48 and 32^ x 64 are generated using DDHMC 
algorithm. The 16^ x 32 configurations were generated using the HMC algorithm except for K = 
0.15775 where DDHMC was also used with different block sizes. The number of thermalized 
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Figure 2: Topological charge for different configurations at volume 32^ x 64, j3 = 5.8 and K = 0.15475 with 
(a) no smearing (b) 20 levels of HYP smearing. 




Figure 3: The Monte Carlo trajectory history for topological charge with unimproved Wilson fermion and 
gauge action for (a) j3 =5.6, volume = 24^ x 48 with a gap of 24 trajectories between two consecutive 
measurements and (b) )3 =5.8, volume = 32^ x 64 with a gap of 32 trajectories between two consecutive 
measurements. 



configurations ranges from 2000 to 12000 and tlie number of measured configurations ranges from 
70 to 500. For topological ciiarge density, we use the lattice approximation developed for SU (2) 
by DeGrand, Hasenfratz and Kovacs [8], modified for SU{3) by Hasenfratz and Meter [9] and 
implemented in the MILC code [10]. To suppress the ultraviolet lattice artifacts, smearing of 
link fields is employed. The comparison of the effect of APE [11] and HYP smearing [12] on 
topological charge is presented in Fig. 1. It is clear that HYP smearing is more effective than 
APE. We used 20 HYP smearing steps with optimized smearing coefficients a = 0.75, a2 = 0.6 
and 053 = 0.3 [12]. In Fig. 2 we show that the topological charges are not close to integers before 
smearing but they are very close to integers after 20 steps of HYP smearing. Fig. 3 shows Monte 
Carlo trajectory history of topological charge at ff = 0.1575,0.15825, j3 = 5.6, volume = 24'' x 48 
and at fc = 0.15455,0.15475, /3 = 5.8, volume = 32^ x 64. There is some sign of trapping only 
at k: = 0.15475 at j3 = 5.8. In table 1 we present, the simulation parameters and the measured 
integrated autocorrelation times for topological susceptibility (T^^j) and topological charge (t^,) 
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Table 1: Lattice parameters and simulation statistics, T is molecular dynamics trajectory length. 



(Two points shown in the presentation were preliminary and are currently going through more 
configuration generation. On the other hand, a few points on the lowest volume have been added, 
especially with different block size of DDHMC). There is some indication that T^^, and T^, decrease 
with increasing K. The T^j at volume = 16^ x 32, j8 = 5.6, K = 0.15775 decreases with increasing 
block size. The T^, at volume = 32^^ x 64, j3 = 5.6 is lower than 24^ x 48 and 16^ x 32 DDHMC 
runs because of the higher block size which results in higher active to total link ratio [13]. In Fig. 
4 we present a few plots for T^^, and T^it determinations. 

Fig. 5 displays four histograms of topological charge distributions, for two values of jS and dif- 
ferent volumes. The topological charge data were put in several bins and the bin widths were chosen 
to be unity centered around the integer values of the topological charges for all the cases. From 
theoretical considerations the distribution of the topological charge is expected to be a Gaussian 
[14]. Since our configurations are reasonably large in number but finite, an incomplete spanning of 
the topological sectors may occur and {Q) may not be zero. Hence we define the susceptibility to 

bez = M(e')-(e)')- 
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Figure 4: The integrated autocorrelation lengths for topological charge and susceptibility at (a) j3 =5.6, 
K = 0.1575, volume = 16^ x 32 (b) j3 = 5.6, k = 0.15775, volume = 24^ x 48 (c) j3 = 5.6, k = 0.15775, 
volume = 32^^ x 64 (d) p = 5.6, K = 0.15815, volume = 32^ x 64 (e) j3 = 5.8, K = 0.15455, volume = 
32^ X 64 (f) j3 = 5.8, k = 0.15475, volume = 32^ x 64. The absolute scales for both the axes are obtained 
by multiplying by 25 in (a), by 24 in (b) and by 32 in the rest. 
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Figure 5: The topological charge distribution for (a) j3 = 5.6, K = 0.1575, volume = 16^ x 32 (b) j3 = 5.6, 
K = 0.1575, volume = 24^ x 48 (c) /3 = 5.8, K = 0.1543, volume = 32^ x 64 (d) j5 ^ 5.8, k = 0.15455, 
volume = 32^ x 64. 



5 



0.2 



f SESAM-T/L (r,, at chiral limit) 

• L16T32 
■ L24T48 

♦ L32T64 



0.15- 



0.1- 



|3=5.6 



2 4 6 8 10 

2 2 

m r„ 



0.002 

0.0015 

l5 0.001 
X 

0.0005 



. L32T64; p = 5.6 
. L32T64; p = 5.8 



0.02 0.04 0.06 0.08 

m (GeV) 



Figure 6: Left Figure: Topological susceptibility versus in the units of ro (at chiral limit) for j3 = 5.6 and 
at lattice volumes 16^ x 32, 24^ x 48, and 32^ x 64 compai-ed with the results of SESAM-TjL collaborations 
[1]. Right Figure: Topological susceptibility versus in the physical units for j3 = 5.6 and j3 = 5.8 at 
lattice volume 32 x 64. The leading order chiral perturbation theory prediction, % — \^ "^g where E is the 
chiral condensate, is also shown for the range 230 MeV <T.t < 290 MeV. 



3. Results 



Fig. 6 (left) shows our results for topological susceptibility versus m^, in the units of Som- 
mer parameter (ro) at the chiral limit, for /3 = 5.6 and at lattice volumes 16^ x 32, 24^ x 48, and 
32^ X 64 compared with the results of SESAM-T^L collaborations [1]. We note that SESAM-T;t:L 
results were presented in [1] by scaling topological susceptibility and m\ by appropriate powers 
of quark mass dependent ro/a. Since r^/a significantly increases with decreasing quark mass, the 
suppression of topological susceptibility is concealed in such a plot. Using the numbers given in 
[1], we have replotted it after scaling by the value of ro quoted at the physical point. The Fig. 
6 (left) clearly shows the suppression of susceptibiUty with decreasing quark mass in the earlier 
SESAM-T;^L data with unimproved Wilson fermion. Our results carried out at larger volume and 
smaller quark masses unambigously establish the suppression of topological susceptibility with de- 
creasing quark mass in accordance with the chiral Ward identity and chiral perturbation theory. In 
Fig. 6 (right) we show topological susceptibility versus nonperturbatively renormalized [15] quark 
mass {rriq) in MS scheme [16] at 2 GeV in physical units for j8 = 5.6 and j3 = 5.8 at lattice volume 
32 X 64. The leading order chiral perturbation theory prediction, % = ^^niq where £ is the chiral 
condensate, is also shown for the range 230 MeV < < 290 MeV. 

In summary, we have addressed a long standing problem regarding topology in lattice simu- 
lations of QCD with unimproved Wilson fermions. Calculations are presented for two degenerate 
flavours. We have presented integrated autocorrelation time for both topological charge and topo- 
logical susceptibility for the two j8 values studied. The effects of quark mass, lattice volume and 
the lattice spacing on the spanning of different topological sectors are presented. The suppression 
of the topological susceptibility with respect to decreasing quark mass, expected from chiral Ward 
identity and chiral perturbation theory is observed. 
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